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i-Sh ■ Abstract 

Q-i 

The goal of this paper is to develop a systematic method of locating the Mueller 
matrices within the class of the matrices of the real group SX(4,R). The main idea is to 
construct the general transformation of the group SL(4, R) (whose real matrices have unit 
determinant) is straightforward, but to analyze the adequacy of such a transformation for 
describing Mueller matrices is highly nontrivial. However, using the technique of Dirac 
matrices, we can quite easy and explicitly describe all the 16 one-parametric subgroups, 
from which, using the all possible products, emerges the whole group SX(4,R). As a 
matter of fact, for these separate 1-parametric subgroups the question of their adequacy 
of describing Mueller matrices becomes sufficiently simple and thus we obtain in each 
case a definite answer. 

1 Elementary 1-parametric generators of the group C7L(4, K) 

It is well known the cornerstone role played in polarization optics by the Mueller matrices 
[3]. Distinguished subsets of the set of Mueller matrices generate group structures, which are 
isomorphic to the group of rotations or to the Lorentz group. 

Since the Mueller matrices are real, of order 4x4, and acting on the 4-dimensional 
real Stokes vector, to investigate the sets of all possible Mueller matrices, one can use the 
parametrization of 4-dimensional matrices which is obtained on the ground of using the basis 
of Dirac matrices, developed in the works [H O [2]. The main goal of this paper is to develop 
a systematic method of locating the Mueller matrices within the class of the matrices of the 
real group SL(4, R). The main idea is the following. To construct the general transformation 
of the group <SX(4, R) (whose real matrices have unit determinant) is straightforward, but 
to analyze the adequacy of such a transformation for describing Mueller matrices is highly 
nontrivial (practically impossible). However, using the technique of Dirac matrices, we can 
quite easy and explicitly describe all the 16 one-parametric subgroups, from which, using 
the all possible products, emerges the whole group SX(4, M). As a matter of fact, for these 
separate 1-parametric subgroups the question of their adequacy of describing Mueller matrices 
becomes sufficiently simple and thus we have reasons to expect in each case a definite answer. 
The explicit form of the Dirac 16-dimensional basis (using the Weyl spinor representation) 
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All these 15 matrices (let us note them as A&) are of Gell-Mann type: - this means they 
have null trace and they are Hermitian too; moreover, the square of each of them is equal to 
I: 

SpA = 0, (A) 2 = /, (A)+ = A, A e {A k : k = 1, 15 }. 
Due to the identity (A) 2 = /, the exponential of each from these matrices A is 

U = e iaA = cos a + i sin ak, det e iaA = +1, a G M. 
We will use the following notations [HE] for some specially chosen six generators A$: 
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Note that we have the commutation relations ctjPk = Pk a j- Due to this commuta- 
tion one can construct nine Abelian 2-parametric subgroups. For instance, e ! " iai e !i ' 1 ^ 1 = 

All the possible products of matrices (jl.2p provide us nine Lie algebra generators: 
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We further specify the explicit form of the 15 elementary 1-parametric unitary transfor- 
mations: 
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and hence we easily get the 15 real 1-parametric 4 x 4-transformations from the group 
SL(4,R): 
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Note that to the generator \q = I corresponds the finite element 
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2 On the (pseudo) Euclidean rotations of Stokes 4-vectors 



A real 4 x 4-matrix M may be considered as being of Mueller type and acting on polarized 
light (completely or partially) if M a i,Sa 



S', 
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M 30 5 + M 3 i5i + M 32 S 2 + M 33 S 3 , 
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In more detailed form, these inequalities look as 



M 00 S + M oiSi + M 02 S 2 + M 03 S 3 > 0, 

(MqoSo + MqiSi + M 02 S 2 + M 03 S 3 ) 2 - (M 10 S + M n Si + M l2 S 2 + M V3 S 3 ) 2 - (2.3) 
(M 20 S + M21S1 + M 22 S 2 + M 23 S 3 ) 2 - (M 30 S + M31S1 + M 32 S 2 + M 33 S 3 ) 2 > 0. 

We recall that So is the intensity / of the light beam, and Si = SoPi = Ipi, where pi is the 
polarization vector. Correspondingly, the above inequalities read as 



Pi + p\ + P% < 1, M oo + Afoipi + M 02 P2 + M 03 P3 > 0, 

(M 00 + Mom + M 02 p 2 + M 03P3 ) 2 - (M10 + M im + M 12 p 2 + M 13 p 3 ) 2 - 

(M 20 + M 2lPl + M 22 p 2 + M 23P3 ) 2 - (M 30 + M 3lPl + M 32 p 2 + M 33 p 3 ) 2 > 0. 



(2.4) 



The problem we face here is very complicated due to the big number of the independent 
parameters - elements of the matrix M a ^. So, we may expect various solutions. In the first 
place, we are interested in Mueller matrix sets which exhibit a group structure. 

The most evident and known such sets are the 3-parametric group of the Euclidean 3- 
rotations, and the 6-parametric group of the pseudo-Euclidean rotations, which form the 
Lorentz group. We shall describe them below. Moreover, we shall consider all the 1- 
parametric elementary generators of the real linear group S , L(4,M). Examination of other 
subgroups is a subject for further concern. 

Let us consider, by using the notation introduced above, the Lorentzian rotations of the 
Stokes 4-vectors. To this end let us specify the following two subgroups: 
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Because the matrices of these two subgroups commute one with each other, we can mul- 
tiply them and, in such a way, we can obtain a new subgroup. Moreover, this new subgroup 
allows us to impose the following constraints to the parameters: 
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These constraints lead to 
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The matrices R(k,k*) determine transformations in the 4-dimensional space with one real 
and three imaginary coordinates. Transition to all four real coordinates is achieved via 
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For real parameters we get Euclidean rotations 
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For complex parameters we get pseudo-Euclidean transformations: 
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1 Note that for the need of polarization optics the transformations having non-unit determinant may also 
be of interest. 
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Let us write below the explicit form of the factors R a (ko, kj) and Rp(kQ, —kj), after perform- 
ing the above similarity tramsformation: 
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According to the afore mentioned considerations, any arbitrary transformation of the 
Lorentz group may be factorized into two mutually commuting and conjugate terms as 

L = K K* = K* K. 

In Physics literature, that fact was firstly noted by Einstein and Mayer (in 1932-1933), while 
constructing the theory of semi-spinors (SKUEj; further, a systematic approach of the Lorentz 
group was developed by Fedorov [2]. 

Let us write below the explicit form of the 2-parametric commuting factors: 
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3 The action of the pseudo-Euclidean rotations on the par- 
tially 

polarized light 

Let us consider the action of pure Lorentz transformations on 4-vectors in the context of 
Stokes formalism, and firstly let us restrict to partially polarized light. The corresponding 
Stokes 4-vectors are analogous to the velocity 4-vectors of massive particles: 

S a = (I,Ip), U a = (U ,U°V); 

Moreover, these 4-vectors behave in the same way with respect to Lorentzian transformations 

cosh /3 — e sinh (3 



L 



-e sinh ft [5ij + (cosh ft — l)ejej] 



1. 
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They act on the Stokes vectors as follows: 



or 



1=1 (cosh (3 — sinh (3 ep), 
I'p' = I [ -sinh P e + p + (cosh /3 - 1) e(ep) ], 



/' = /(cosh p — sinh pep), 
^ - sinh P e + p + (cosh p — 1) e(ep) 
cosh /3 — sinh /3 ep 

Let us specify them in the following several special cases: 
1: 

n t> t u « / ~ Sinh ^ e + P 
ep = 0, /=/cosh/3, p= — — , 

cosn p 

where the intensity of the light beam increases and the polarization degree decreases. 
2: 

p = +p e, 1=1 (cosh p — sinh P p), 

where the polarization vector changes according to the law 

, , — sinh P + cosh P p 
^ ^ cosh P — sinh P p 

Therefore the polarization degree varies as 

, 1 — p" 1 tanh p 
^ 1 — tanh P p 

Note that when p > 0, this decreases, and when P < 0, then it increases. In particular, there 
exists the "rest reference frame" where the partially polarized light becomes natural light: 

L = L (p ,n), tanh Po = P, 

p' = 0, I' = I (cosh po - sinh tanh p ) = — ^— , (3.2) 

cosh po 

and here the polarization degree equals to zero. 
3: 



-p e, I' = /(cosh P + sinh /3p), 



p' = p > c _ sinh + cosh g c (33) 
cosh /3 + sinh P p 



The polarization degree varies according to the law 

1 + tanh p 
1 + tanh /3 p ' 



where for P > increases, while for /3 < it decreases. Again, there exist an analogue of the 
"rest reference frame" Lq, in which the polarization degree equals to zero. 
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Let us further consider the role of the relativistic ellipsoid in optics. We start with the 
simplest case, e = (0,0, 1), while the formulas (3.2) give 



I' = I (cosh (3 — £>3 sinh j3), 
Pi 



, cosh /3p3 — sinh (3 
^ 3 cosh P — y>3 sinh j3 ' 
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Due to the basic property of the Lorentz matrices, the following identity holds: 

r'2/i '2\ j2(-\ 2\ i '2 1 — p 2 

/ (1-p ) = /(l-p) or 1-p = (cosh/g _ p3sinh/g)2 . 

Therefore, the polarization degree transforms according to 

1-p 2 
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cosh /? + sinh f3 p 3 



cosh /3 + sinh /3 p 3 
and (13. 6h reduces to the form 

p' 2 = 1 - (1 - p 2 )(cosh /3 + sinh f3 p' 3 ) 2 . (3.7) 
Let us show that this equation (|3.7p describes an ellipsoid. Indeed, (|3.7h can be re-written as 
Pi 2 + Pa 2 + P3 + (1 - P 2 ) 2 cosh /3 sinh /3 p' 3 + (1 - p 2 ) sinh 2 /3 p 3 2 = 1 - (1 - p 2 ) cosh 2 /3 , 



or 



+ P2 2 + (cosh 2 /3 - p 2 sinh 2 /3) 
= p 2 cosh 2 /? — sinh 2 /3 + 



, (1 — p 2 ) sinh (3 cosh /3 
^ 3 cosh 2 (3 — p 2 sinh 2 /3 
(1 -p 2 ) 2 sinh 2 /3 cosh 2 ^ 

cosh 2 (3 — p 2 sinh 2 /3 



and finally we obtain the ellipsoid equation 

Pi + p'i + (cosh 2 (3 - p 2 sinh 2 /3) (p 3 + 7) 2 

where 



cosh 2 (3 — p 2 sinh 2 /3 ' 



(3.8) 



1 



(1 — p 2 ) sinh /3 cosh /3 
cosh 2 f3 — p 2 sinh 2 /3 



cosh 2 (3 -p 2 sinh 2 /3 = cosh 2 /3(1 - p 2 ) + p 2 > 0. 



Thus, the surface which has the form of a sphere (p 2 = p 2 ) will transform under the 
action of a Mueller matrix of Lorentz type into ellipsoid (|3.8p . 
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This result can be extended for more general Mueller matrices of Lorentz type of arbitrary 
orientation: 

, . , p - e sinh (3 + (cosh g - 1) e (ep) 

1=1 (cosh j3 - smh (3 (ep) , p = — —— . 

cosh p — smh pep 

Considering the identity 

/WWd-/) - i-f-^z-jlZ^—p, (3.9) 

then, excluding the variable p, we get: 

p' + e sinh (3 + (cosh j3 — 1) e(e p') 



P 



cosh j3 + sinh /3 e p' 
or 

cosh P - sinh /3(ep) = — — — -. 

cosh p + smh p ep' 

From (|3.9p we get that 

1 — p 2 = (1 — p 2 ) (cosh /3 + sinh (3 ep' 2 . 
This describes an ellipsoid equation with the orientation governed by the vector e. 

4 The action of a Lorentzian transformation on completely 
polarized light 

Now let us consider the action of Lorentzian transformations on completely polarized light 
(the analogous of the isotropic 4- vectors in Special Relativity): 

S a = (IJn), n 2 = l. (4.1) 

With respect to transformations of kind 



L = 

the Stokes 4-vector transforms as follows 



cosh f3 — e sinh f3 

e sinh (3 [5ij + (cosh (3 — l)eiej] 



, . , -sinh ^e + n + (cosh /3 - 1) e(en) 
1=1 (cosh p — smh p en), n = 



cosh f3 — sinh (3 en 
Let us specify now several special cases: 
1: 

n t' t x, a / n - sinrl P e 
e n = , / = / ch p, n = — . 

cosh p 

2: 

e - + n n ' _ ~ sinh (3 n + n + (cosh (3 - 1) n _ ^ 
' cosh f3 — sinh f3 
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3: 

ji+13 i smn n + n + (cosh [3 — 1) n ^ (4 2) 

' cosh P + sinh /3 

The known non-existence of the "rest reference frame" for massless particles means, in the 
context of polarization optics, that one cannot completely transform the polarized light into 
a natural one. 

Note that under the action of Mueller transformations of Lorentz type, the degree of the 
polarization of light does not change, that is p' = p = 1. 



5 On the deformations of the Stokes 4-vectors 

Now, let us consider in detail the Mueller matrices of diagonal form: - these describe simple 
deformations of the Stokes 4-vectors. Firstly, let it be the variant Uq(i\): 



M = U (i\) 

















e~ x 














e~ x 














e~ x 



s„ 



e- x S a . 



The inequalities (|2.ip hold good: 



S n 



e~ x S > 0, 



5" 
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e~ 2X Sl > 0. 



The action of such Mueller matrices on the light is given by the relations 

I'- X I, 

We further consider the variant ?7^(iA): 



Pi =Pi 



M 



e X S , 



e x 













e~ x 














e x 














e~ x 


s[ = 


e~ x 


Si, 


S 



e X S 2 , 



e~ X S 3 . 



The restrictions (|2.ip take the form (the first inequality always holds good): 



e A S > 0, 



S 



e 2X {Sl -Si)- e~' ZA {S'f + S'£)> 0. 



The second one is equivalent to 



o2 , o2 
^4A \ °1 + °3 



> 



c2 



o2- 
^2 



Since the initial light obeys the restriction 

c2 _ q 



sf - s| > 



1 > 



s 2 + si 



C2 _ o2 
°0 '-'l 



2 ' 



(5.1) 
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where the inequality (2.2b) will be true when A is positive, that is 

A G [0,+oo) 

The action of these Mueller matrices on the light is given by the relations 

I' X I, p' 1 = e~ 2X Pl , p' 2 = p 2 , 
Now, consider the variant U^(i\): 



P's 



e" 2A P3. 



M 



S' 
S' 2 



e~ x 














e" A 














c x 














e x 


s[ = 




Si, 


S' 2 



S' 3 = e x S 3 , 



e~ 2X (S 2 - S 2 ) - e 2X (S 2 + S 2 ) > 0. 
The last inequality is equivalent to 

o2 q2 



> e 



A G (-oo,0]. 



sf + Sf 

The intensity of the light and the degree of polarization transform as follows: 
I'- X = I, p[=Pi, 
Now, consider the variant U% (iX): 



P2 = e 2X P2, 



P 3 = e 2X p 3 , 



M 

S'o 
S' 2 



Ef (iA) - 

-2\/q2 



e~ x 













c x 













e x 










e~ A 




s[ = 


~-e x 


Si, S' 2 


= e 


sl)- 


e 2X 


(S 2 + S 2 2 )>0. 



S' 3 = e" A ft, 



The last inequality gives 



c2 q2 
°0 °3 > g 4A 



A e (-oo,0]. 



Pi =e 2A Pi, 



P2 = e2X P2, P 3 = IP3 



(5.2) 



(5.3) 



+ SI 

The intensity of the light and the degree of polarization transform as follows: 
l' x = I, 

Instead of the four diagonal Mueller transformations described above, one may introduce the 
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following simpler four elementary deformations: 



E-2 



E, 



\ 

C A 
















1 










n 
U 


U 


1 
1 


n 
U 


J 


() 


o 


o 


1 




1 
















e A 
















1 





J 











1 




e~ 


X 
















e A 
















e A 
















e" A 


1 
















1 
















1 





J 











e A 





^' = I, Pi =Pl, P2 =P2, P 3 



P3, 



I' = I, Pi = eX Pi, P'l = P2, p' 3 = P3, 



I' = I, P'l = Pi, P2 = eX P2, P3 = P3, 



I' = I, Pi = Pi, p'i = P2, p' 3 = e X P3- 



(5.4) 



It is obvious that the above four matrices determine a 4-parametric Abelian subgroup of 
the mutually commuting transformations. 



6 On other subgroups of 5L(4, R) 

Among the generators (1.5) one can emphasize the following triplets: 



and 





K 


= {Ai 


= OL\Pl 


B 2 


= a 2 /3 2 ,C 3 


= a 3 (3 3 }, 




L 


= {d 


= OLlPz, 


A 2 


= a 2 pi,B 3 


= a 3 f3 2 }, 




M 


= {Bi 


= aif3 2 


c 2 


= a 2 /3 3 ,A 3 


= a 3 Pi}, 


K' 


= { 


-C x = - 


-ai/3 3 , - 


-B 2 


= -a 2 f3 2 , - 


-C 3 = -a 3 (3 3 } 


V 


= { 


-J3i = - 


-uifo, - 


-A 2 


= -a 2 Pi, ~ 


-B 3 = -a 3 fi 2 } 


M' 


= { 


-Ai = - 


-otiPx, - 


-c 2 


= -a 2 /3 3 , - 


-B 3 = -a 3 f3 2 } 



(6.1) 



(6.2) 

For all 3- vectors of the above generators, we can accordingly construct multiplication laws of 
the same type: 

rir 2 = -r 3 , r 2 r! = -r 3 , r\r 2 - r 2 r! = o, 

and similar ones provided by cyclic permutations. Within each triple, the generators commute 
one with each other. This means that on the base of each triplet of generators one can 
construct 1-, 2-, 3-parametric Abelian subgroups. 
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Let us mention that using the 15 generators ctj, fy, Ai,Bi,Ci, one can build 20 isomorphic 
to SU(2) Lie algebras: 



(ai, 0:2,03), 
(ax,A 2 ,A 3 ), 
(ai,B 2 ,B 3 ), 
(oi,C 2 ,C 3 ), 
(Pi,Bi,Ci), 
(^i,&,Ci), 
(Ai,J3i,&), 



(A,/3 2 ,/3 3 ), 

(Ai,a 2 , A 3 ), 

( J Bi,a 2 ,-B 3 ), 

(Ci,a 2 ,C 3 ), 

(Pi,B 2 , C 2 ), 

(A 2 ,f3 2 ,C 2 ), 

(A 2 ,B 2 ,fo), 



(A 1 ,A 2 ,a 3 ), 

(B 1 ,B 2 ,a 3 ), 

(C 1 ,C 2 ,a 3 ), 

(Px,B 3 ,C 3 ), 

(A 3 ,fa,C 3 ), 

(A 3 ,B 3 ,p 3 ). 



(6.3) 



All these say that the reserve of the possible subgroups of SX(4, R) is rather large, and that we 
may expect to separate some of them as being appropriate to define Mueller-type subgroups. 

7 Examination of 16 elementary one-parametric subgroups in 

SX(4,R) 

Let us consider the remained 12 one-parametric subgroups in 5X(4,R): 
(7.1) 

Variant L7f (</>): 



M = U?(<j>) 



or with the use of the block form 



COS (f> 


sin (ft 








— sin</> 


COS (j) 














COS <fi 


— sin 








sin 6 


cos q 



M 



k + k 3 ki + k 2 

h -k 2 fe - ^3 

^0 + ^3 ^i + Z 2 

^1 — ^2 ^0 — ^3 



n + n 3 
ni - n 2 
m + m 3 
mi — m 2 



ko = cos <f> 
mo = cos </> 
n = 
/o 







fcl 
mi 

m 
/1 










fc 2 

m 2 



ni + ra 2 
n - n 3 
mi + m 2 
m - m 3 

= + sin (j) , 
= — sin 4> , 
"-2 = 0, 

/ 2 = o, 



The restrictions (2.1) give 



So > 



c2 
°0 



fc 3 
m 3 

"3 
h 



c-2 
D 2 



S, 2 > 



cos 20(5o 



S?) + sin20 25 5i 



(cos ^>5o + sin 4>S\) 2 

°2 — D 3 



> 
> 



cos 4>So + sin 05i > , 

(— sin 4>Sq + cos 4>S±) 2 + S 2 + S : 

. 



(7.1) 
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With respect to the variables 

tan 4> = x , pi = a , p 2 = b , p 3 = c, 
the above inequalities take the form 

a sin cj) + cos > , 
^(l-a 2 ) + -^2a-6 2 -c 2 > 0. (7.2) 

1 + X z 1 + X z 

The sets of solutions of the first inequality can be illustrated by the following formulas and 
by Fig. 1 (the interval <j> £ [l,27r] is divided into four parts): 

1 

a > , cos(£ > 0, <fr G {/, IV}, x>--\ 

a 

a > , cos0 < 0, € {II, III}, x < ; 

a 

a < , cos0 > 0, G {I, IV}, x< — ; 

a 

a < , cos0 < 0, € {/J, III}, x>—. (7.3) 

a 

The quadratic inequality (j7.2|) is equivalent to 

(a 2 - 1 - b 2 - c 2 )x 2 + 4ox + (1 - a 2 - b 2 - c 2 )) > 0. 
whose solution is x G [xi,X2], where x\, x 2 are the roots of the equation 



2a 1 - a 2 - b 2 - c 2 



r 2 — 2r — = 

6 2 + c 2 + l-a 2 6 2 + c 2 + l-a 2 ~ 



These are the real numbers 



2a - y^a 2 + (1 - b 2 - c 2 - a 2 ) (b 2 + c 2 + 1 -a 2 ) 
Xl " 6 2 + c 2 + 1 - a 2 ' 

2a + y^a 2 + (1 - b 2 - c 2 - a 2 ) (b 2 + c 2 + l^ 2 ) 
X2 ~ b 2 + c 2 + 1 - a 2 ■ ( } 

In the case of completely polarized light, the formulas become much simpler: 

6 2 + c 2 + a 2 = l, xi2 = jo, 



1-a 2 



and two different possibilities should be distinguished 



2a 
l~d 2 

2a 

~a l 

a = 0, x\ = X2 = (degenerate case) 



a > 0, 0<x<- ^ {^ in h 



a < 0, ^<x<0 {II, IF}; 
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x = tan (j) 



a > x = tan (ft 



a < 



7 II III IV 



I II III IV 



Figure 1: Solutions of inequalities (7.1.4) 
For a > 0, the joint solution of both inequalities in (|7.2[) are possible only in the quadrant 

I: 

2a 

a > 0, < x < 



1-a 2 ' 

For a < 0, the joint solution of both inequalities in (|7.2p is possible only in the quadrant 

IV: 

2a 

a < 0, = < x < 0. 

1 — cr 

Let us investigate now the more complicated case of partially polarized light. Since for 
any value of a the root x\ is negative, and the root X2 is positive, turning to Fig. 1, we see 
that the joint solutions of both inequalities in (|7.2p are possible. 



(7.2) 

Variant Ug((j)) 



M = (-0) = 



COS (j) 





s'mcf) 








COS 4> 





sin(/!> 


— sin (j) 





COS 4> 








— sin0 





COS (j) 



or in block form 



M = 



k + h 


h 


+ k 2 


n + 


»3 


ni + n 2 




ki - k 2 




-h 


ni - 


»2 




- n 3 




k + h 


h 


+ h 


m + m 3 


mi + 7712 




h-h 


lo 


-h 


mi — 


m 2 


m 


- ni 3 




ko = cos 




kl 


= o, 




= 0, 


h-- 


= o, 


niQ = cos 




mi 


= 0, 


?n 2 


= 0, 


m 3 -- 


= 0, 


no = sin 


0, 


n\ 


= 0, 


»2 


= 0, 


n 3 -- 


= 0, 


Iq = — sin 




h 


= 0, 


/2 


= o, 


h- 


= 0. 



(7.5) 
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The restrictions (|2.ip lead to 

S >0 



C2 



si 



s 2 - s 2 > 0, 



Ql'2 



s' 2 > s' 2 2 + s[ 2 



3 • 



5q = cos ^5*0 + sin 05*2 > 0, 
(cos (pSo + sin 4>S 2 ) 2 — (cos0Si + sin 0<S 3 ) 2 — 
(— s'm(j)So + cos 4>S2) 2 — (— sin^Si + cos 4>Ss) 2 > =^ 
cos20(S 2 - S 2 2 ) + sin 20 25 5 2 - Sf - Sg > . 
In the variables a,b,c,x, these inequalities take the form 

< 1, cos 6 + b sin > 0, 



1 



:(l-6 2 ) + 



2;r 



a 2 + 6 2 + c 2 
26 - a 2 - c 2 



0. 



(7.6) 

1 + x z 1 + x z 

They differ from the previous ones from (|7.2[) only in formal notation, and therefore the 
results will be much the same. 

(7.3) 

Variant Ug(<p) 



cos 








sin0 





COS 


— sin0 








sin 


COS 





— sin0 








COS 



or in block form 



M 



k + k 3 
h - k 2 
k + h 
h-h 



ki + k 2 
k - h 
h + h 
h - h 



ko = cos 0, 
mo = cos 0, 
n = 0, 
Z = 0, 



h 
mi 

m 
h 



n + n 3 
ni - n 2 
m + m 3 
m i — m 2 
= 0, 



0. 
0. 
0, 



ni + n 2 
n - n 3 
mi + m 2 
m - m 3 

k 2 = 0, 
m 2 = 0, 
n 2 = sin 0, 
Z 2 = sin 0, 



h 



0. 
0. 
0, 
0. 



(7.7) 



The restrictions (|2.ip lead to 



S > 0, 



c2 



5? 



o2 
D 2 



S 2 >0, 



S'v 

sin 0S 3 )' 



o 



> 
> 



Si 



(cos 0So 

=► cos 20(5^ - 5|) + sin 20 2S S 3 
In the variables a, 6, c, x, the inequalities take the form 

cos + c sin 



cos 0So + sin 05 3 > , 
(— sin05o + cos 0<S 3 ) 2 + S 
S\ 



S\ 



1 



1 + x 2 



(1 



c 2 ) + 



2.r 



1 + x 2 



2c -a" 



> 
> 



0. 



(7.8) 
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These differ from the previous ones from (|7.2p only in formal notation, and therefore the 
results will be much the same. 



(7.4) 



Variant U 



M 



COS(j) 


sin <j) 








— sin</> 


cos 4> 














cos 4> 


sin (ft 








— sincf) 


COS (j) 



or in block form 



M 



ko + k 3 k\ + k 2 n + n 3 n\ + n 2 
k\ - k 2 k - k 3 m -n 2 n - n 3 
lo + h h + h m + m 3 mi + m 2 
h ~h h~ h mi -m 2 mo — m 3 
k$ = cos</>, k\ = 0, k 2 = + sin</>, 
mo = cos<^>, mi = 0, m 2 = + sin<^>, 
n = 0, m = 0, n 2 = 0, 

l = 0, h = 0, h = 0, 



^3 

h 



0. 
0. 
0. 
0. 



The restrictions (|2.ip give 



So 
S'o 

(cos </>Sq + sin (pSi) 2 



> 



S, 



o 



sf > s 2 + s 2 , 



cos 0So + sin 4>S\ > , 

(— sin^So + cos^Si) 2 + S 2 + S : 

} , C-2 



> 



> s 2 + s 2 



cos 2(j){S^ - S() + sin 2<£ 2S Si 
With respect to the variables a, b, c, x, the inequalities become 

cos d> + a sin <A > 



1 + x 2 



2x 



1 + x 2 



2a - b 2 - c 2 > 0. 



(7.9) 



(7.10) 



These differs from the previous ones from (17. 2j) only in formal notation, and therefore the 
results will be much the same. 



(7.5) 

Variant ui( 



COS (j) 





sin (j) 








cos 4> 





— sin 


— sin</> 





COS (j) 








sin <h 





COS 
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or in block form 



M 



k + k 3 ki+ k 2 n + n 3 n x + n 2 
ki - k 2 k — k 3 ni - n 2 n Q — n 3 
+ ^3 h + h mo + m 3 mi + m 2 
h —h lo - h mi - m 2 m - m 3 
ko = cos (j), k\ = 0, 
niQ = cos (j), mi = 0, 
no = 0, n\ = 0, 
l = 0, h = 0, h = 0, h 



k 2 
m 2 
n 2 
h 



0, 
0, 
0, 
0, 



m 3 



0, 
0. 



"3 



The restrictions (|2,ip give 



sv. 



Ql'2 

b o o 



'2 



5i 
S'o 



> 
> 



o2 
°3 > 



5l 2 



sin i 
sin i 



S 2 



(7.11) 



Si > 0, 



(sin^>5i + cos 4>S 3 ) 2 
coe20(Sg - Sf) + sin 20 2S 5 2 - S 2 - Sf 



> 
> 



(— sin <j)So + cos 05"2 

In the variables a, b, c, x the above inequalities take the form 

< 1, 



cos 4>So + sin (pS 2 > 0, 
(cos (f)So + sin 4>S 2 ) 2 — (cos 0Si 

0. 



sin^Ss)' 



a 2 + b 2 + c 2 



COS ( 



6 sin 6 > 0, 



1 



1 + 



:(1 



6 2 ) + 



2;r 



1 + x 2 



26 



> 



They differ from the previous ones from 
will be much the same. 



0. (7.12) 
2]) only in formal notation, therefore the results 



(7.6) 

Variant U 3 (4>) 



M = Ul{^) = 



COS (j) 








sin 





cos 4> 


sin (j) 








— sin0 


cos 4> 





— sine 


b 





cos 



or in block form 



M = 



k + k 3 


ki + k 2 


n + n 3 


ni + n 2 








ki - k 2 


k - k 3 


n\ - n 2 


n - n 3 








h + h 


h + h 


mo + m 3 


mi + m 2 




1 




h-h 


h — h 


m\ — m 2 


mo - m 3 








ko = cos 




h = o, 


k 2 = 0, 




h 


= o, 


mo = cos 




mi = 0, 


m 2 = 0, 


m 3 


= 0, 


n = 


o, m 


= sin (j), 


n 2 = 0, 




™3 


= 0, 


k = 


o, h = 


— sin (fi, 


i 2 = o, 




h 


= 0. 



(7.13) 
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The restrictions (|2.ip give 



So >0, 



c2 
^0 



5? 



5 2 - Si > 0, 



5'' 2 



o 



S' 2 > S' 2 2 + S' 2 



3 • 



S"q = cos 05q + sin 053 > 0, 



) 2 - (cos05i + sin05 2 ) 2 - 
- (— sin05o + cos 4>S 3 ) 2 > 



(cos 05o + sin 
(— sin 4>S\ + cos 05 2 ) 2 - 

cos 20(5 2 - 5|) + sin 20 25 5 3 - 5f - 5 2 2 > 0. 
In the variables a, b, c, x the inequalities take the form 

< 1 



1 



1 + x- 



■(l-c 2 ) + 



2x 



l + x 2 



a 2 + b 2 + c 2 



2c- 



cos + c sin > 0, 



> 0. 



(7.14) 



These differ from the previous ones from (|7.2|) only in formal notation, therefore the results 
will be much the same. 



(7.7) 

Variant U^(-if3) 



U 2 A l 



cosh j3 



sinh j3 





cosh j3 
— sinh f3 






— sinh f3 
cosh /3 




sinh (3 



cosh /3 



or in block form 



M 



k - 
ki- 
lo- 

h- 

k 

mo 



h 
h 



k\ + k 2 
k - k 3 
h + h 
h — h 
cosh /3, 
cosh /3, 
n = 0, 
Z = 0, 



no + n 3 
n\ - n 2 
mo + m 3 
mi — m<i 

h = o, 

mi = 0, 

m = o 
h = o 



»2 
h 



n\ + n 2 

n - n 3 
m\ + m 2 
mo - m 3 

k 2 = 0, 
m 2 = 0, 
= — sinh /3, 
= — sinh /?, 



m 3 
h 



0. 
0. 
0. 
0. 



e^(5o + 53) + e-^(5o-5 3 ) > 0, 



The restrictions (|2.ip give 

cosh /35o + sinh /3S 3 > =>■ 
which holds good for any (3. At the same time we have 

(cosh /35 + sinh ,55s) 2 - (cosh /3Si - sinh /35 2 ) 2 - 
-(- sinh pSi + cosh /35 2 ) 2 - (sinh (3S + cosh f3S 3 ) 2 > 
cosh 2/3(5 2 + 5f) + 2 sinh 2f3S 1 S 2 > 0. 



(7.15) 



c-2 
°0 



Si 
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In the variables a, b, c and y = tanh /3, y £ (—1, +1), the second inequality takes the form 

1 - c 2 - ^4(« 2 + b 2 ) + j^2ab > 0, 



or 

-y 2 (a 2 + 6 2 + 1 - c 2 ) + 4a6y + (1 - a 2 - 6 2 - c 2 ) > 0. 

The coefficient of y 2 is positive, so the inequality holds good in the interval between the roots 
of the corresponding quadratic equation, that is 

y e [2/1,2/2], y\ < tanh p < u 2 , 



2ab - y/4aW + (1 - a 2 - b 2 - c 2 )(a 2 + b 2 + 1 - c 2 ) 
2/1 ~ a 2 + b 2 + 1 - c 2 <0 ' 

2ab + v / 4a 2 6 2 + (1 - a 2 - b 2 - c 2 )(a 2 + b 2 + 1 - c 2 ) /w . 

2/2 = 2 , , 2 , -i 2 > °- 7 ' 16 

a z + o z + 1 — c z 

Let us check the necessary condition y 2 < +1; we have 

y / 4a 2 6 2 + (1 - a 2 - b 2 - c 2 )(o 2 + b 2 + 1 - c 2 ) < (a 2 + & 2 + 1 - c 2 ) - 2a6, 

which after being squared, leads to 

4a 2 6 2 + (1 - a 2 - 6 2 - c 2 )(a 2 + 6 2 + 1 - c 2 ) < 
< (a 2 + b 2 + 1 - c 2 ) 2 + 4a 2 6 2 - 4a6(a 2 + b 2 + 1 - c 2 ). 

This is equivalent to 

< 2(a 2 + b 2 ) - 4ab < (a - 6) 2 . 

Therefore, the relation y 2 < +1 is satisfied. 

Now let us verify the second necessary condition y\ > —1: we have 



- ^4a 2 6 2 - (1 - a 2 - b 2 - c 2 )(a 2 + b 2 + 1 - c 2 ) > -(a 2 + 6 2 + 1 - c 2 ) - 2a6, 
which yields 

4a 2 6 2 + (1 _ a 2 _ fe 2 _ c 2 )(a 2 + 6 2 + x _ c 2) < 

< (a 2 + 6 2 + 1 - c 2 ) 2 + 4a 2 6 2 - 4a6(a 2 + b 2 + 1 - c 2 ). 
This is equivalent to 

< 2(a 2 + b 2 ) - 4ab <^ < (a - ft) 2 , 

and hence the relation yi > — 1 holds good. 

For the case of completely polarized initial light, we have a 2 + 6 2 + c 2 = 1 and the formulas 
become much simpler: 

y G [yi,2/2], 2/1 < tanh /3 < y 2 , 



ab — V a 2 b 2 ab + V a?b 2 

yi = a 2 + b 2 ' 2/2 = a 2 + b 2 ' 
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so that 



ab > 0, Ui = 0, 



ab < 0, 



V2 



2ab 



< +1; 



2ab 



> -1, 



2/2 = 0. 



(7.17) 



(7.8) 

Variant E#(i/S) 



cosh /3 — sinh /3 

cosh (3 — sinh /? 

— sinh /3 cosh (3 

— sinh (3 cosh (3 



or in block form 



M 



k + k 3 ki+ k 2 n + n 3 ni+ n 2 

k\ - k 2 k - k 3 ni -n 2 n - n 3 

^0 + ^3 h + h mo +rn 3 mi + m 2 

h ~h h~ h m\ -m 2 mo — to 3 

ko = cosh f3, k\ =0, k 2 = 0, k 3 = 0, 

too = cosh /3, toi = 0, to-2 = 0, to 3 = 0, 

no = — sinh /3, m = 0, n 2 = 0, n 3 = 0, 

Z = - sinh /3, Z x = 0, l 2 = 0, Z 3 = 0. 

The restrictions (|2.ip lead to 

cosh /3So — sinh /3S2 > holds good for all /3, 
(cosh PS Q - sinh /3S 2 ) 2 - (cosh /3Si - sinh j3S 3 ) 2 - 
-(- sinh PS + cosh /3S 2 ) 2 - (- sinh /35i + cosh /3S 3 ) 2 > =► 
S o ~ s i ~ cosh 2 P( S i + 5 i) + 2 sinh 2/3 1 S 1 S 3 > 0. 
In the variables a, b, c the second inequality becomes 

1 - b 2 - cosh 2(3(a 2 + (?) + sinh 2/3 2ac > 0. 
This differs from (I7.2p only by notations, so the results will be much the same. 

(7.9) 

Variant Uf(i/3) 



cosh j3 

cosh {3 

— sinh /3 

— sinh (3 



— sinh j3 
— sinh f3 
cosh /3 

cosh j3 



(7.18) 
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or in block form 



M 



k + k 3 
k\ - k 2 

k + h 

h-h 

ko = cosh (3, 
mo = cosh (3, 
n = 0, 
l = 0, 



k\ + k 2 
ko - k 3 
h + k 

^0 — ^3 



n + n 3 
n\ - n 2 



n\ + n 2 
n - n 3 



m + m 3 mi + m 2 
mi — m 2 mo — m 3 



h 



fci = 0, 


k 2 


= 0, 


^3 


= 


mi = 0, 


m 2 


= 0, 


m 3 


= 


- sinh /3, 


»2 


= 0, 


"3 


= 


- sinh /3, 


/2 


= 0, 


/3 


= 



The restrictions (|2.ip give 



cosh (3 So — sinh /35 3 > 0, 
which holds good for all f3. At the same time we have 

(cosh (3S - sinh /3S 3 ) 2 - (cosh /3Si - sinh /35 2 ) 2 - 
-(- sinh f3Si + cosh /3S 2 ) 2 - (- sinh /3S + cosh (3S 3 ) 2 > 



c-2 



c-2 
^3 



cosh 2/3(S 2 + Si) + 2 sinh 2/3SiS 2 > 0. 
In the variables a, b, c the second inequality becomes 



1 - c 2 - cosh 2/3(a 2 + b 2 ) + sinh 2/3 2a6 > 0. 
This differs from (I7.2p only in notation, so the results will be much the same. 

(7.10) 

Variant U^ifi) 



U 3 B (i(3) 



cosh (3 
sinh j3 







sinh f3 
cosh (3 










cosh f3 
— sinh (3 






— sinh f3 
cosh /3 



or in block form 



M 



ko + k 3 ki + k 2 
h - k 2 k - k 3 
lo + h h + h 
h — h h — h 
ko = cosh f3, 
mo = cosh (3, m\ 
n = 0, 
/ = 0, 



n + n 3 ni + n 2 

ni -n 2 no- n 3 

mo + m 3 mi + m 2 

mi — m 2 mo — m 3 

k\ = sinh (3, k 2 = 0, 

- sinh f3, m 2 = 0, 

ni =0, n 2 = 0, 

h =0, h = 0, 



k 3 
m 3 
n 3 
h 



0, 
0. 
0. 
0. 



The restrictions (|2.ip give 



(7.19) 



(7.20) 



cosh (3So + sinh j3S\ > 0, 
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which holds good for all f3. At the same time we have 



(cosh /3S + sinh pSi) 2 - (sinh /3S + cosh /3Si) 2 - 
-(cosh (3S 2 - sinh /3S 3 ) 2 - (- sinh /3S 2 + cosh (3S 3 ) 2 > =► 
S$ - S 2 - cosh 2/3(S 2 + Sf) + 2 sinh 2/3S 2 S 3 > 0. 

In the variables a, b, c the second inequality becomes 

1 - a 2 - cosh 2/3(6 2 + c 2 ) + sinh 2/3 26c > 0. 

This differs from (|7.2p only in notation, so the results will be much the same. 

(7.11) 



Variant Uf(i/3) 



cosh /3 


sinh /3 






cosh /3 


— sinh j3 



sinh /3 


cosh j3 






— sinh j3 


cosh j3 



or in block form 



M 



k + k 3 ki + k 2 n + n 3 n\ + ra 2 
fcl - & 2 ^0 — &3 - «2 - ™3 

+ ^3 h + h m + m 3 mi + m 2 
^i — — ^3 mi - m 2 m - m 3 

&o = cosh /3, fei = 0, k 2 = 0, 
mo = cosh /3, mi = 0, m 2 = 0, 
n = 0, n\ = 0, n 2 = 0, 
l = 0, Zi = 0, l 2 = 0, 



"3 
Z 3 



A: 3 = 0, 
m 3 = 0, 
= sinh /?, 
= sinh /3. 



The restrictions (|2,ip give 

cosh /3S + sinh /3S 2 > 0, 
which holds good for all /3. We also have 

(cosh (3S + sinh f3S 2 ) 2 - (cosh /3Si - sinh /3S 3 ) 2 - 
-(sinh /3S + cosh (3S 2 ) 2 - (- sinh ps x + cosh /3S 3 ) 2 > 
=> S$ - Si - cosh 2/3(5? + Si) + 2 sinh 2/3SiS 3 > 0. 

In the variables a, 6, c the second inequality becomes 



eP(l + 6) + e - P(l-6) > 0^ 
!> 2 - cosh 2/3(o 2 + c 2 ) + sinh 2/3 2ac > 0. 



(7.21) 



(7.22) 



This differs from (|7.2p only by notation, so the results will be much the same. 
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(7.12) 



Variant Ug(-0) 



-%P) 



cosh P 
— sinh /3 





— sinh f3 
cosh /3 










cosh (5 
— sinh (3 






— sinh (5 
cosh f3 



or in block form 



M 



k + k 3 ki+ k 2 

ki - k 2 k - k 3 

lo + h h + h 

h — h h — h 

ko = cosh /3, ki 
mo = cosh /3, mi 
n = 0, 
l = 0, 



n + n 3 
ni - n 2 



n 



n>2 



mo + m 3 mi + m 2 

mi — 777-2 m — m 3 



sinh /3, 
sinh /3, 

77,1 = 0, 
il =0, 



A: 2 

777 2 

"-2 
h 



0. 
0. 
0. 
0. 



k 3 = 0, 
m 3 = 0, 
77 3 = 0, 

h = o. 



The restrictions (|2.1|) give 

cosh /3So — sinh f3S\ > , 
which holds good for all f3. We also have 

(cosh (3 So - sinh /3Si) 2 - (- sinh /3S + cosh /3Si) 2 - 
-(cosh /3S 2 - sinh (3S 3 ) 2 - (- sinh /35 2 + cosh /35 3 ) 2 > 
S 2 - S 2 - cosh 2/5(5| + 5|) + 2 sinh 2/3S 2 S 3 > 0. 

In the variables a, b, c the second inequality becomes 



1 



cosh 2(3(b 2 + (?) + sinh 2/3 26c > 0. 



It differs from (|7.2p only in notation, so the results will be much the same. 



(7.23) 



8 Varying the degree of polarization of the light 

Let us examine how the degree of polarization of the light changes for the above 12 special 
one-parametric Mueller transformations. The corresponding characteristics may be given by 
the difference 



(a '2 + b '2 + c '2 } 



(a 2 + b 2 + c 2 ) = D. 
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The D-entities for the first six cases (7.1)— (7.6) are 



(7 1 1 


n 

LJ — 


(a- 


x) 2 + {b 2 + c 2 )(l + x 2 ) 


- a 2 


- b 2 


_ r 2 




(1 + ax) 2 




(7 0\ 
{1. 1) 




(b- 


x) 2 + {a 2 + c 2 ){l + x 2 ) 


— a 


h 2 
— 






(1 + bx) 2 


— C 


(7.3) 


D = 


(c- 


x) 2 + (a 2 + b 2 )(l + x 2 ) 


-a 2 


-b 2 


-c 2 




(1 + ac) 2 




(7.4) 


D = 


(a- 


x) 2 + {b 2 + c 2 )(l + x 2 ) 


-a 2 


-b 2 


-c 2 




{l + ax) 2 




(7.5) 


D = 


(b- 


x) 2 + {a 2 + c 2 ){l + x 2 ) 


-a 2 


-b 2 


-c 2 




(1 + bx) 2 




(7.6) 


D = 


(c- 


x) 2 + {a 2 + b 2 ){l + x 2 ) 


-a 2 


-b 2 


-c 2 




(1 + ac) 2 





(8.1) 

It will be sufficient to consider in detail only one case, e.g., (7.1). In this case, we get: 

(a-x) 2 {b 2 + c 2 ){x 2 + 1) 2 l2 2 
(l + xa) z (l + xo) z 

The degree of polarization decreased when 

D< 0, (a - x) 2 + (6 2 + c 2 )(x 2 + 1) < p 2 (l + xa) 2 , 

or when 

x 2 (l + b 2 + c 2 ) - 2ax < x 2 a 2 p 2 + 2xap 2 , 

and this yields 

D<0, x 2 (l - a 2 ) - 2ax < 0. 
Thus, the light is depolarized if 

2a 

D < 0, < x < 7T , a > 0; 

1 — GT 

-D < 0, — ^ < x < , a < 0. (8.2) 
1 — a z 

In contrast, the degree of polarization of the light increases if 

x 2 (l - a 2 ) - lax > 0, (8.3) 

that is, when 

2a 

D > 0, x > = , a > 0; 

1 — a z 

D > 0, x < — % , a < 0. (8.4) 
1 — a z 

The light does not change the degree of polarization if 

2a 

D = => x = - ? . 

1 — a 2 
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Figure 2: The line D(a,x) = 

This fact is illustrated by Fig. 2. 

Now let us examine the other six cases (7.7)— (7.12). 
The case (7.7) (we use the notation tanh f3 = y) gives 

_ (a cosh j3 — b sinh /3) 2 + (-a sinh f3 + b cosh /3) 2 + (sinh j3 + c cosh /3) 2 

(cosh /3 + c sinh /3) 2 



_ a 2 _ ft2 _ c 2 = y^-^) 2 + (6-ay) 2 + (c + j/) 2 _ fl2 _ 6 2 _ ^ 



(8.5) 



(a - 6y) 2 + (6 -ay) 2 + (c + j/)' 
(1 + cy) 2 

The degree of polarization decreases when 

D < 0, (a - ft?/) 2 + (6 - ay) 2 + (c + yf < (a 2 + o 2 + c 2 )(l + cy) 2 , 

or when 

y 2 (a 2 + b 2 + 1) + 2y(c - 2a6) < y 2 c 2 p 2 + 2ycp 2 , 

and this leads us to 

D < , y 2 (l - c 2 )(l + p 2 ) - 2y[2ab - c(l - p 2 )] < 0. 
The solutions have the form 

D < 0, 2ab - c(l - p 2 ) > 0, < y < _ 2 - — [2a6 - c(l - p 2 )] , 

D < 0, 2a6-c(l-p 2 ) < 0, 2 [2ao - c(l - p 2 )] < y < 0. (8.6) 

(1 — c z J(l + p z J 

The degree of polarization increases when 

D > , (a - 6y) 2 + (6 - ay) 2 + (c + y) 2 > (a 2 + b 2 + c 2 )(l + cy) 2 . 



J We use the constraint a 2 + b 2 + c 2 — p 2 



29 



This gives 

y 2 (l - c 2 )(l + p 2 ) - 2y[2ab - c(l - p 2 )] > 0. 

The solutions are 

D > 0, 2ab-c(l-p 2 ) >0, y>- 2 [2afe - c(l - p 2 )] , 

(1 — c z J(l + p z J 

D > 0, 2ab-c(l-p 2 ) <0, y< (1 _ c2) 2 (1+p2) [2a6-c(l-p 2 )] . (8.7) 

The light does not change the degree of polarization if 

D = 0, y= (1 _ c2) ( 1+p2) [2^-c(1-p 2 )]. 

The formulas become considerably simpler when the initial light is completely polarized. 
All the six variants of this type are characterized by the following D-entities: 

(7.7) D = (a -^ )2+ / b - a f + ( ^ + c)2 -a 2 -6 2 -c 2 , 

(1 + cy) 2 

(7.9) D = ^-by? + (h-ay) 2 + (y~c) 2 _ a 2_ b 2_ c 2^ 

(1 - cy) 2 

(7.8) D = {a ~ CV? + ( ( i C jg 2 2 + {V ~ 6)2 - a 2 - b 2 - c 2 , 

(7.11) D = (Q -^ )2 + (c -^ )2 + ( ^ + 6)2 -a 2 -6 2 - C 2 , 

(1 + 6y) 2 

(7.10) D = (fe -^ 2+ / C -^ 2 + ^ + a)2 -a 2 -6 2 -c 2 , 

(1 + ay) 2 

(7.12) D = (b-cy) 2 + (c-by) 2 + (y-a) 2 _ q2 _ # _ ^ (g 8) 

(1 - ay)^ 

It is obvious that the differences between these relations consist of notations only, so the 
derived results will be similar. 

Examination of other subgroups (see Section 7) is a subject for further concern. 
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